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Abstract 

We prove an identity on Hermitian random matrix models with 
external source relating the high rank cases to the rank 1 cases. This 
identity was proved and used in a previous paper of ours to study the 
asymptotics of the top eigenvalues. In this paper, we give an alter- 
native, more conceptual proof of this identity based on a connection 
between the Hermitian matrix models with external source and the 
discrete KP hierarchy. This connection is obtained using the vertex 
operator method of Adler and van Moerbeke. The desired identity 
then follows from the Fay-like identity of the discrete KP r vector. 

1 Introduction 

The subject of this paper is an identity between a Hermitian random matrix 
model with external source of rank m and m such models with external 
source of rank 1. This identity allows us to reduce the asymptotic study of 
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"spiked Hermitian random matrix models" of rank higher than 1 to that of 
the models of rank 1. This reduction formula was used in |2] to evaluate the 
limiting fluctuations of the top eigenvalue(s) of spiked models of arbitrary 
fixed rank with a general class of potentials. In [5] we gave a direct proof 
of this identity using the formula of |3] on the Christoffel-Darboux kernel. 
Here we give a different, more conceptual proof using the relation between 
the random matrix model with external source and discrete KP hierarchy. 
We show how the general results of Adler and van Moerbeke [1] on the 
construction of solutions of discrete KP hierarchy can be used on the partition 
functions of the Hermitian matrix model with external source. 

We now introduce the model. Let W{x) be a piecewise-continuous func- 
tion on M. We assume that W{x) is nonnegative, has infinite support, and 
vanishes sufficiently fast as — )■ oo so that converges. Let A he a. d x d 
Hermitian matrix with eigenvalues ai, ■ ■ ■ , a^. We call A the external source 
matrix and W{x) the weight function. Consider the following measure on 
the set Hd of d x d Hermitian matrix M: 

P{M)dM := — ^ det{W{M))e^'^^^'''^dM, (1) 

Zd{ai, ■ ■ ■ ,ad) 

where W{M) is defined in terms of a functional calculus and 

Zd{a^, ■■■ ,ad):= [ det{W{M))e^''^^^UM (2) 

is the partition function. Note that the partition function depends only on the 
eigenvalues of A, but not its eigenvectors. It is also symmetric in ai, ■ ■ ■ , a^. 
The Harish-Chandra-Itzykson-Zuber integral formula [6], [8] implies that 

Zd{ai, ■■■ ,ad) 
Cd 



n 



-/ det[e^^'''Y _^ n i^!^-^^)IlW{X,)dX, (3) 

- «i) ' i<i<k<d j=i 



for a constant Cd which depends only on d. There is a similar formula for 
the density function of the eigenvalues. 
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If some of eigenvalues of A are zeros, we use the following short-hand 
notations: For m < 



Zd{ai, ■ ■ ■ , am, 0, ■ ■ ■ , ) 




(4) 



The main theorem of this paper is about the following expectations. For 
-E C M, s G C, and m < d, we define 



(£rf(ai, • • • ,am;E; s) := E 



- sxE{\j))PiM)dM, 



(5) 



where Ai, ■ ■ ■ , are eigenvalues of M and the expectation is with respect 
to the measure ([1]) when the eigenvalues of the external source matrix A are 
oi, ■ ■ ■ ; (^m, 0, ■ ■ ■ ,0. We note that the last integral is the partition function 



d—m 



with new weight (1 — sxe{x))W{x) divided by the original partition function 
with the weight W{x). This observation will be used in the proof later. When 
s = 1 the above expectation is a gap probability. Set 

(£d(ai, 



<£d(ai, 



1 Om, S , 



,am;E;s) 



(6) 



^d{E;s) 

Let Pj{x) be orthonormal polynomials with respect to W{x)dx. For a real 
number a, define the constant 



Pj{s)e'''W{s)ds. 



(7) 



The main result is 
Theorem 1.1. We have 

€d(ai, • • • , ttm', E; s] 



det [Td-j{ak)^d~j+iiak; E; s)]^^^^ 



det[r,_,(a,)].%i 



(8) 



if ai, - ■ ■ ,am are distinct and non-zero. 
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Since both sides of (jS]) are analytic in a^, the above identity still holds 
when some of aj are the same or equal to zero if we interpret the right-hand 
side using THopital's rule. 

The term ^a-j+iidk', E; s) on the right-hand side of ([8]) is given by ([6]) 
with the rank m = 1 and the only non-zero eigenvalue of the external source 
A being Cfc, and the dimension is changed to d — j + 1. Hence the identity 
relates the rank m case to m rank 1 cases. 

In [1], the asymptotics ^d-j+i{cik', E; s), the rank 1 cases, were obtained. 
The quantities Td-j{ak) were also analyzed asymptotically in the same paper 
as an intermediate step. In the subsequent paper [S] the higher rank cases 
were then analyzed asymptotically using the identity ([H]). The main result 
was that when the potential, assuming that it is real analytic, is convex to 
the right of the right-end point of the support of the equilibrium measure, the 
phase transition behavior of the fluctuations of the top eigenvalues is same as 
in the Gaussian unitary ensemble. Otherwise, new types of jump transitions 
are possible to occur. A characterization of possible new transitions was also 
obtained. 
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2 Proof 

Set A(x) := ni<j<fc<d(^fc - Xj) = det{x']~^)'^f,^^ for x = {xi, ■ ■ ■ ,Xd). Recall 
that, setting a = (ai, ■ ■ ■ , a^) and A = (Ai, ■ ■ ■ , A^), 

det ^ s4\)s4a) 

where the sum is over all partitions k, with at most d parts. Here denotes 
the Schur polynomial and denotes the number of parts of partition k. 
This identity can be proved as follows. First, from the Andreief 's identity [2] 
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(equivalently the Cauchy-Binnet formula), 



det [eM.];^,. det E ^ E -^J-^" . (10) 



>- n=0 



1 - det [a;^] det [Aj'j 



Since the summand is symmetric in n^'s and equals to zero when two of the 

oo 

summation indices are the same, we can replace ^ by 

ni,--- ,nii=0 0<7id< - <r!.2<ni 

Now set Kg := Uq — d + q. Then the summation condition becomes Ki > ■ ■ ■ > 
> 0, i.e all partitions with at most d parts. The identity ([9]) follows by 

recalling the classical definition of the Schur polynomial ^^(a) = ^p^^ . 

Inserting (|9]) into (|3]), we obtain the Schur polynomial expansion of the 
partition function: 

where 

« d 

G^:= s,{X)A{X)'l[w{X,)dX,. (12) 

Using the classical definition of the Schur function, the determinantal form 
of A (A), and the Andreief's identity [2], we obtain 

G^ = d\- det [M^^+d-p+q.,Y Mj := / x^W{x)dx. (13) 

Here Mj are the moments of the measure W(x)dx. 

We insert (fT3!l into (fTTIl and use the Jacobi-Trudi identity, 3^(0) = 
det [/iKp-p+g(a)]p 5=1, where /ij(a) denotes the complete symmetric function. 
The sum is over the partitions k = (ki, ■ ■ ■ , k^) where > ■ ■ ■ > > 0. 
We set jp := Hp + d — p. Then ji > ■ ■ ■ > jd > 0. Since the summand is 
symmetric in jp's and vanishes when two indices are the same, we arrive at 
the formula 

, ^ det [M,^+.-i]',,=idet [/^,,-d+,(a)],^,,=i 
^d(2i,--- ,ad) = C'd 2^ . (14) 
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In the below, the notation t = (ti, ^2? ■ ■ ■ ) denotes a sequence of variables. 

2 3 

We also use the notation [c] = (c, y, y, ■ ■ ■ ) for the evaluation of t by powers 
of c. The notation [ci] + [C2] + ■ ■ ■ + [c^] stands for the evaluation of t obtained 
by substituting tj = c^/j- 

Following ini Definition 6.1], we define the so-called "elementary Schur 
polynomials" hj{t) by the generating function 

00 

J2 /ij(tK = eSr-i*^"'\ (15) 

j=-oo 

If t = [fli] + [02] + • • ■ + [ad] (i.e. tj = Yli=i'^ilj)^ hj{t) is the complete 
symmetric function in ai, ■ ■ ■ , a^, which we denoted earlier by hj{a). This 
abuse of notations is unfortunate but in the below we only use the definition 
h{t) given in (fT5l) . 

Now define the formal power series in t = (ti, ^2, ■ ■ ■ ) 

Z,(t)-=^ y det[M,^+,_i]^^,^,det[/i,^_,+,(t)]^,^,, ^^^^ 

where Cd '■= -^j^- This definition is equivalent to the formula (26) of [12]. 
Then 

Zd{ai, ■■■ ,ad) = Zd{[ai] H h M)- (17) 

Setting some of the parameters to be zero, we also have for m < d 

Zd{ai, ■ ■ ■ , am) = Zd{[ai] H h [a^]). (18) 

We now show that Zd{t) solves the discrete KP hierarchy following the 
vertex operator construction of the general solutions due to Adler and van 
Moerbeke. We then show that a general property of the vertex operator 
solution implies an identity of which the identity (j8]) is a special case. Note 
the definition (fT6|) does not require the assumption that the weight function 
W{x) is nonnegative. In the next two subsections we drop this assumption. 

2.1 Discrete KP r vector 

We observe that: 
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Proposition 2.1. Let Z^it) be defined in fll6p and set Zq^I) := 1. Then the 
sequence 

{■■■ , Z^i-t), Zi(-t), Zo(t), Zi(t), Z2{t), ■■■) (19) 

constitutes a discrete KP r vector. 

Discrete KP r vectors are solutions to a system of differential-difference 
equations in the discrete KP-hierarchy. Adler and van Moerbeke established 
several characterizations of discrete KP r vectors in [1] . Here we use two of 
them: the vertex operator characterization (see Proposition 12.21 below) and 
the Hirota bilinear identity characterization (see (155]) below). 

Remark 2.1. Each components of a discrete KP r vector is a KP r function. 
The fact that Zii{t) is a KP r function for each G N was proved in [12]. In 
[7] it was proved further that for each d G N, ^d(t) is a so-called 1-KP-Toda 
r function and moreover these 1-KP-Toda r functions are derived from the 
same Grassmannian structure. It should also be possible to prove the above 
proposition from this fact. 

The vertex operator is a differential operator defined by (see [H Formula 
(0.22)]) 

5^4.Mexp -5^^— . (20) 

fc=i ^ ^ fc=i 
The vertex operator acts on a formal power series /(t) as 

X(t,z)/(t) = eS^-*^^V(t-r^])= ( £ /ifc(t)z'=y(t-[z-i]). (21) 

Adler and van Moerbeke found a very general construction of a discrete KP 
r vector from one KP r function and a sequence of measures using the vertex 
operator: 

Proposition 2.2 ([H Theorem 0.3]). Let r(t) 6e a KP r function. Let {■ ■ ■ , 
v^i{z)dz, VQ{z)dz, Vi{z)dz, ■■■) 6e a sequence of arbitrary measures. Then 
the infinite sequence (■ ■ ■ , T_i(t), ro(t), ri(t), ■ ■ ■ ) defined as ro(t) := r(t) anc? 

rrf(t) := X(t, z)ud^^{z)d^ ■ ■ ■ Q X{t, z)iyo{z)dz^ r(t), (22) 
r„,(t) := X(-t, z)iy_d{z)dz^ ■ ■ ■ (^j X(-t, 2)z/_i(2)rf^^ r(t), (23) 
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for c? G N, forms a discrete KP r vector. 



Remark 2.2. The original statement of [H Theorem 0.3] assumes that the 
measures Vk{z)dz are defined on M. However, it is easy to check that the 
proof to Proposition 12.21 in [1] holds almost verbatim if we change M into C 
(or more restrictively the unit circle {zgC||2;| = 1} that we are going to 
use in this section). 

Proof of Proposition \2. li We set r(t) = ^o(t) := 1. This is trivially a KP 
r function. Hence Proposition 12.11 is proved if we can construct a sequence 
of measures (■ ■ ■ , h'-i{z)dz, h'o{z)dz, vi{z)dz, ■ ■ ■) such that Zdit) equals the 
right-hand side of ([22]) (resp. ([23D) with r(t) = 1 for > (resp. d < 0). 
Define the measure on the circle {zGC:|2;| = l}as 

^^^z)dz := Izl)!^ y Mi±l,d-j-i^^^ d = 0, 1, 2, ■ ■ ■ , (24) 
27riCd+i ^ j! 

where Cq := 1. We also define 

Udiz) := u^d-iiz) rf=-l,-2,---. (25) 

We now show that these measures satisfy the desired property. 

For d = 0, since Zo{t) = 1,^ implies that X{t, z)Zo{t) = J2T=~oc> hkit)z^- 
Hence we find from a direct evaluation of the integral using the Cauchy in- 
tegral formula (the integral is over the unit circle) that 

<fxit,z)Zo{t)uo{z)dz ^ 1 = z^{t) (26) 



from the definition (11611 . 

We now consider d > 0. From fl2T]) . 

X{t,z)Ut) = -^( £ h,{t)z') 
d\Cd \u / 



g det[M,^+,-i]^,^^,det[/i,^-,+,(t - [z~^])]i,=. 



(27) 
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Note that from we have Ej°l-oo^j(t " = e^^^^^*^"^^"'' - 

(l ~ ■7) J2'jL-oo hj{t)wK Comparing the coefficients of , we find that 

Substituting ( 125]) into (|27|) . we can derive, after some algebra, that 



... ^ _ — n ; n -J 



M. 



30 



M. 



n 



,jd=0 k=0 

ji—d+1 (t) 



jd+d-l 



h 



■jd-d. 



't) h, 



jd-d+l 



hjdiK 



(29) 



Note that there is one more summation index jo and the determinants are 
of (i + 1 hj d + 1 matrices. Then from (1291) and (ITB]) . and noting that the 
variable z appears only in the last column of the ffist matrix in the (l29l) . we 
can check directly using the Cauchy integral formula that 



X{t, z)Zd{t)vd{z)dz = Z,+i(t), > 0. 



Successive applications of the relation (130!) imply that for all > 0, 



X(t, z)vd-i{z)dz 



which is same as (|22|) . 

Finally, ^ and ([SID imply that 



X{t, z)vQ{z)dz I Zo(t 



(30) 



(31) 



Zd{-t) = (^j X{-t,z)p.d{z)d^---(^j X{-t,z)u^^{z)d^ Zo(t). (32) 



This is same as (1231) . Hence the proposition is proved. 



□ 



2.2 Fay-like identity 

An importance property of discrete KP r vector is that its components sat- 
isfy a Hirota bilinear identity (see lU Theorem 0.2(iii)]). (Adler and van 



9 



Moerbeke, moreover, showed that the Hirota bihnear identity actually char- 
acterizes the discrete KP r vector.) For the discrete KP r vector fll9p in our 
situation, this identity becomes 

-i. / Z,,{i - [z~'])Z,,+,{t + [^"i])eSr^^fe-*.>^-^'^-^-id^ = (33) 

for all di > (^2 > 0. Here the formal integral of a formal Laurent series is 
defined by 

^£^( E^«^-^)^- = «- (34) 

We now show that this Hirota bilinear identity implies a Fay-like iden- 
tity fHOj) . Such a derivation of a Fay- like identity from the Hirota bilinear 
identity was obtained in the Toda lattice hierarchy by [11] and we adapt this 
approach. 

We take the special choices di = d, d2 = d — 2 and t = t + [a] + [b] in fl5^ . 
Then the factor e^j"^!^*^"*-')''''^;'^!"^^-! gq^^jg 2gE^i(x+T)^\ After re-writing 
the sum in the exponent as — log(l — az) — log(l — bz), and using the simple 
identity 7- — "ff, , , = rAr(r-^ -^) + i, we find that 

{l—az)(l—bz) {a—b)z^l—az 1—bz' z' 



00 



z^^^T^T>' = —J-— y a^z^ - —-1— y Vz^ + ^. (35) 
aia - b)z ^ b(a - b)z ^ abz 



Using this, ( 1331) implies that 

a / _ , 1 . / , ^ A dz b / _ , , , / v-^ \ dz 



-L«--)(E-)^-iL«-')(E 



a^z^ 



2tii V ^ J z 27ri J^^^ V ^ / ^ 

a — b f ^/ ,dz 
Q{z 



(36) 

where Q is defined by 

Q(ti;) := Z,{t + [a] + [b] - [w])Z,_,{t + M). (37) 

Observe that the Laurent series of Q{w) consists only of non-negative powers 
of w. Hence 

00 

Q{w) = J2qnW'' (38) 

n=0 
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for some qo,qi, - ■ ■ ■ Thus, from the integral 

dz 

= 2^qna" = Q{a) 
j=0 ^ " n=0 (39) 

= Zrf(t + [a])Zrf„i(t + [6]). 



Similar evaluations of the other integrals of (139|) imply the following Fay-like 
identity: 



aZd(t + [a])Z,_i(t + \h\) - bZdit + [6])Z,_i(t + [a]) 
= {a-b)Zd{t + [a\ + [b])Zd-i{t). 



(40) 



In the remaining part of this section, using identity ( HOj) we prove that: 



Proposition 2.3. For any d > m > 1 and ai 

Zdit + [ai] + ■ ■ ■ + [a^]) 1 



ZAt) 



det 



lai, ■ ■ ■ , Or 



Zrf+i_j(t + [ofc]) 



i,fc=i 
(41) 



where Am(ai, ■ ■ ■ , ctm) := ni<j</fc<m(% ~ '^fc)- 

Proof. After dividing the identity (1401) by Zd(t)Zrf_i(t), we obtain 



^d(t 



[^]) 



^d(t) 



det 



Zd(t) " Zd(t) 
id,i(t+[a]) irf„i(t+[fe]) 



2d-l(t) 



L(t) 



(42) 



This is the identity ( 1411) when m = 2 for all (i > 2. 

We now prove the general case using an induction in m. Suppose that fl41l) 
holds with m < m — 1 for all (i > m — 1 and ai, ■ ■ ■ , Om-i G C. We are to 
prove that it holds with m = m for all > m and ai, ■ ■ ■ , G C. For this 
purpose, we set a = ai, b = a^, and t i— >■ t + [02] H — ■ + [o-m-i] in (|42|) . After 
pulling out the denominators of the entries of the determinant outside, we 
obtain 



(ai - am)Zd{t + [ai] H h [am])Zd_i(t + [02] H h [am,-i]) 



det 



aiZdit + [ai] 
Zrf_i(t + [ai] 



[a^-i]) arnZdit + [02] 
[flm-l]) Zd„i(t+[a2]- 



[a-m]) 

[flm]) 



(43) 
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Let us call the entries of the last determinant Aij, i,j = 1,2. First we 
consider Aij on the first row. From the induction hypothesis, 

-| m— 1 

m-l-jZct+l-jit + [flfc]) 



A 



11 



ai 



Za{t) Am-i{ai, ■ ■ ■ ,a.m-i) 



det 



Zd+i-j{t) 



(44) 



j,k=i 



If we multiply 02 ■ ■ ■ Om-i on both sides and bring the factor ai ■ ■ ■ a^-i inside 
the determinant, we find that 



^2 ■ ■ ■ O-m-l 
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ZJt) A„_i(ai, 



det [5jfc]i<j<„,_i,i<fc<„_i (45) 

[ai, ■ ■ ■ , am-i) — 



where 



m-j ^d+l-j(t + [flfc]) 
Zrf+i_j(t) 



(46) 



Note that the power of is changed to m — j from m — 1 — j of (H4l) . 
Similarly, we find that 

Au 1 



^2 ■ ■ ■ flm-l 



det [Bj,]^ 



<j<m-l,2<k<m 



(47) 



with the same definition of fij^. Note the difference of the indices of the 
determinant from (HSil . 

Now we consider in the second row. The induction hypothesis implies 
that 

m—l 



121 



1 



Zd-i{t) A„_i(ai, ■ ■ ■ , ttm-i) 



det 



a. 



_i_jZd-j{t + [ak]) 



't) 



(48) 



j,fc=i 



Note that d is changed to (i — 1 in the determinant from If we shift the 
index j by j — 1 in the determinant, we can write the above as 

A21 1 



^rf-l(t) Am-l(ai, ■ ■ ■ , Om-l 



■ det [5 



ifeJ2<j<m, l<fc<m-l • 



Similarly, 



A 



22 



Zrf_i(t) Am_i(a2, 



det [5 



(49) 



(50) 
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Consider the matrix B of size m whose entries are Bj^, j,k = 1, ■ ■ ■ , m. 
Let denote the matrix of size m — 1 obtained from B by deleting the row 
a and the column b. Then the determinants in fH5|) . f H7|) . f H9|) . and f l50|) are 
the determinants of the matrices ^B™, Bj^, B^, and respectively. Hence 
we find that 



-l)Am-l(a2, ■ ■ ■ , flm) 



Zrf_i(t)Zrf(t) 
det K] det [Bl] - det [S^] det [S™] . 



det 



■^11 ^12' 
^21 ^22 



(51) 



Now the Desnanot-Jacobi identity (see e.g. [TOl Proposition 10]) implies that 
the above equals det [B] det [-8^'^] where BI'^ is the matrix of size m — 
2 obtained by deleting the rows l,m and the columns l,m from B. The 
determinant det [S] is precisely the determinant in fHT]) with m = m. On 
the other hand, 



det [Bl'-] 

^2 ■ ■ ■ a-m-1 



det 



det 



-i-j Zd+i-j{t 



-1 1 



2'd+i-j(t) 

_2-j^d-j(t + [flfc+l]) 



(t) 



J i,fe=2 
m-2 



i,fc=l 



(52) 



The last determinant is precisely the determinant in fHTl) with m = m — 1, 
d replaced by c? — 1, and the complex numbers given by 02, ■ ■ ■ , am-i- The 
induction hypothesis implies the identity 



det [Bl'-] 



■^d-i(t + [^2] + ■ ■ ■ + [a m-l]) 



Combining (jS]), f lFIjl . and and noting that 

^^-2(^2, ■ ■ • , flm-l) 



(53) 



(ai — am)A.m-i(ai, 



-l)Am-l(a2, 



Am(ai, 



(54) 



we obtain pil) with m = m. Hence the induction step is established and the 
proposition is proved. 

□ 
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2.3 Proof of the theorem 

Now we complete the proof of Theorem ll.li 

Proof of Theorem For any subset E eM. and s G C, consider the new 
weight function We,s{.x) := W{x){l — sxe{x)). Let Z^'^{ai, ■ ■ ■ , a^) be the 
partition function ([2]) with W replaced by We,s- We also use a similar short- 
hand notation as (jl]). Then 

Taking t = (0, 0, ■ ■ ■ ) in fHTl) and recalling ([T7|) . we find 

-j -Z^d+i-i(Qfc) 



Zd{ai, • ■ ■ , a„i) 1 J , 

— ^ V ^ A — 7 V 



m- 



(56) 



Note that this holds for any weight function W. We substitute W i— )■ We,, 
in fl56|) and divide this identity by fl56|) with ly. From this we obtain 



^ , det [a^ ^Zd+i„j(afc)Cd+i_j(afc;E;s), .^^^ 

(Bd{ai,--- ,a^;E;s) = -— ■ . (57) 

det [Ofc 'Zd+i-j{ak)\.j^^^ 

We now consider the terms a^~-' Zfi+i-j{ak)- For any dimension /, if a; = a 
and ai = ■ ■ ■ = ai_i = 0, applying THopital's rule to (E]), we have a formula 
of the partition function Zi{a): 

Zi{a) = ^ , / det [V^] det [Ar^] H W{X,)dX, (58) 

where V = {V^jy..^^, with = \{~^ for j = 1, ■■■,/- 1 and l/^z = e'^^\ Let 
be orthonormal polynomials with respect W{x)dx. By using elementary 
row operations, 

Zi{a) = ! det [V] det [p,-i(A,)] n^(^^)A- (59) 
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where V = (Vij)^,,.^^, with % = Pj-i{\i) for j = I,-- - , / - 1 and Vu = 
6°"^% and C'l is a new constant which depends only on / and W. Using the 
Andreief's formula and the fact that pj are orthonormal polynomials, we 
obtain 

Zi{a) = / e^^pi_,{X)W{X)dX = l\C',a''^'ri_^{a). (60) 
Inserting this into (!57l) . we obtain Theorem 11.11 □ 
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